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We propose a coherent framework allowing to deal with many-body effects in dense ultracold 
Fermi gases in the presence of a Feshbach resonance. We show that the simple effect of Pauli 
exclusion induces a strong modification of the basic scattering properties, leading in particular 
to an energy dependence of the effective scattering length on the scale of the chemical potential. 
This results in a smearing of the Feshbach resonance and provides a natural explanation for recent 
^ experimental findings. 
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• Besides the continuing progress in understanding Bose Einstein condensation in ultracold bosonic atomic gases, the 
T-H ' exploration of similar dense [1] fermionic systems has received recently a strong impetus from experiments reaching 

the strongly degenerate regime with mixtures of fermions in two different hyperfine states [2]. A main purpose in 
exploring these systems is the search for a transition to a BCS superfluid [3]. A particularly fascinating possibility 
. of experimentally controling the strength of the interaction, already demonstrated in Bose systems, has been used in 
c/j ■ these last experiments. It consists in working in the vicinity of a Feshbach resonance [4] where the scattering length, 
^ directly linked to the interaction, depends strongly on the applied magnetic field. R has been emphasized recently 
[5] that the strong interaction obtained in the vicinity of this resonance, together with the Bose condensation of the 
^ \ molecules corresponding to the underlying bound state, could be of major interest for the BCS transition. 

In the present paper we set up a theoretical framework for handling many-body effects in the presence of a Feshbach 
resonance. As a first consequence we show that the simple effect of Pauli exclusion explains qualitatively that the 
Feshbach resonance is strongly smeared for degenerate gases, as it appears experimentally [2]. Here we will restrict 
^ ourselves for simplicity to the case where the (quasi) bound state responsible for the resonance is not thermally 
occupied, which corresponds to a negative effective scattering length. Similarly we will consider only the normal state 
^ ' although our formalism can be extended to the superfluid state. 
<^ \ First we model the Feshbach resonance in the following way. The Feshbach resonance [4] is actually produced by 
. the somewhat complex interplay between the spin (electronic and nuclear) and orbital degrees of freedom. However 
' in order to have an understanding of the physics linked to this resonance, it is convenient to use a simple modeling 
04 where the spin degrees of freedom are no longer involved. One can think that the resonance occurs because there is 
a (quasi) bound state caused by the existence of a deep well in the atomic interaction potential at short distances. 
One will indeed have a bound state, corresponding physically to the formation of a molecule, if this well is essentially 
isolated from the large distance region by a high barrier with very small transmission probability. Clearly this model 
behaves phenomenologically exactly in the same way as a Feshbach resonance [6] . 
"j^ On the other hand we are only interested in the effect of this molecular state on the scattering of two atoms, and 

' we will ignore the effect of the other atoms on these two when they are close together. This is reasonable since these 
I other atoms will be most of the time far away and will not perturbe the two atoms we consider. In other words we 
make use of the fact that the gas is dilute on the scale of the molecular state. Naturally this means that we ignore 
Ch , for example the possibility that three atoms are close together, which would require additional ingredients in the 
O 1 description. 

We further simplify our description by separating the possible distance between the two particles in two domains. 
Either they are far away and interact by the long range part of the potential (this range will be denoted by index 1). 
Or they are quite close and interact through the deep well of the potential (this is noted by index 2). More precisely 
one can define the boundary between the two domains as being at some distance R, large compared to the molecular 
' size but small with respect to interparticle distance. Then instead of describing exactly the transition from domain 
1 to domain 2, we assume that a term in the Hamiltonian gives rise to a matrix element producing this transition. 
We note that this problem is actually quite similar to the one raised by tunneling of electrons between two metallic 
electrodes through an insulating barrier and our approach is at the level of the tunneling Hamiltonian [7] . It is 
actually known that this simple modeling can be improved up to the exact problem, and that many-body effects can 
also be included in this theory [8] . However this does not seem necessary in the present case and our simple modeling 
should be a quite enough. 

To be more specific now consider first the problem with the center of mass of the two atoms at rest, so we deal 
with a one body problem with interparticle separation r and reduced mass to,- = to/2. The parts of the Hamiltonian 
H corresponding respectively to domain 1 and 2 are: 
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where Vi(r) and V2(r) are the interaction potential at long and short distance respectively. For example Vi(r) gives 
rise to a scattering length Ubg while V2(r) produces the molecular state. Then we describe hopping between the two 
domains by the non-local term: 

ifi2(r,r')=i(r,r') (2) 

with H21 = -ffi2- Since the purpose of t{r, r') is just to make the particle cross the boundary between domain 1 and 

2, we can assume it to be short-ranged around this boundary. 

Let now G^xCr, r', w) be the propagator corresponding to Hn, at frequency u> (we take h = 1) and similarly 
G22(r, r', w) the propagator corresponding to i?22- Physically they describe the motion either at long or short distance, 
without the possibility to hop between the two domains. We treat now exactly the effect of the hopping term. We 
introduce the full propagator G(r, r', w) corresponding to H. If r and r' are large (we call again Gn the corresponding 
propagator), the particle can propagate either by staying in domain 1 (this is described by Gji) or, after a stay in 
domain 1, by making a first hop to domain 2 and then propagating by any means back to domain 1. This leads to 
the following equation between operators: 

Gil = Gil + G'iifG2i (3) 

where a same frequency is understood in all the propagators. We have a similar relation for propagation between 
domain 2 and 1, except that hoping is now required (because G21 =0): 



Carrying Eq.(4) into Eq.(3) gives: 



G21 = Gl^t^Gii (4) 



Gn = G°i + G°itG^2t+Gii (5) 



which is an integral equation for Gn. 

This can now be simplified if we take into account that a single bound state of H22 , corresponding to the Feshbach 
resonance we are interested in, is relevant. All the other bound states are supposed to be very far away. This allows 
us to make, in the energy range we are interested in, a single pole approximation for G22: 

GO,(r,r',.) = ^£)|:(£l (6) 

where (p(r) is the wavefunction of the bound state and Eq its energy. This makes Eq.(5) explicitely soluble because 
G22 becomes basically a projector on the bound state. One finds: 

Gn = G?! + —J-^Gl.tl^ X ^|ttG°i (7) 

where SEq =< ip\t^Giit\(f > is a complex quantity. 

Let us assume for simplicity that there is no background scattering, i.e. the long distance potential Vi is zero. 
In this case Gji is just the free particle propagator. The corresponding T-matrix is then given by the last term in 
Eq.(7) without the Gfi operators. For the scattering we are interested in, we look for matrix elements between plane 
waves with very small wavevectors compared to the molecular scale. Since (p(r) and t{r, r') are short-ranged we can 
as well take these wavevectors to be zero. This leads to a numerator equal to |wp with w = J drdr't{r, r')ip{r'). The 
denominator gives a pole for uj = E with E = Eq + SEo, corresponding to the resonance produced by the bound state. 
The real part KeE gives the physical energy of the resonance, which is the one actually measured experimentally. 
So we do not have to worry about calculating Re6Eo. The imaginary part gives the width of the resonance due 
physically to the possible decay, induced by t, of the molecule into two atoms. Introducing the Fourier transform Gj^ 
of the free particle propagator, this imaginary part will come from ImG^ = —itS{u! — ek) with ek = k"^ /2m,r, physically 
linked to the density of final states for the decay. Since we are concerned with low energy lo, the wavcvector must 
be small and the matrix elements coming from t\ip > in the above expression of 5Ei^) can again be evaluated for zero 
wavevector, which introduces again \w\'^ . Finally we obtain for the corresponding scattering amplitude: 

/(k) = \-, (8) 
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to be evaluated on the shell lu = /2mr. In this expression we have set 7 = m,.|wp/27r. We find in particular in 
Eq.(8) that Im/^^(«;) = — k as required by unitarity. Evaluating Eq.(8) at zero energy gives the scattering length 
a = —^/loq. Strictly speaking the Feshbach resonance corresponds to the situation where the above resonance occurs 
at zero energy w = 0. This occurs for wq = 0, i.e. for an infinite scattering length. Experimentally loq is controlled by 
the applied magnetic field. Naturally this result for the scattering amplitude is well known [6] , as well as this general 
way of modeling the Feshbach resonance [4] as a simple switch between molecular state and diffusion states. We have 
just reformulated this approach in a way which lends itself conveniently to generalization in order to include many 
body effects. 

We turn now to the case of a dense Fermi gas and assiimc again for simplicity that there is no background scattering. 
As we have already mentionned we take advantage that the gas is dilute on the molecular scale to neglect the effect 
of the other atoms when two atoms are scattering due to the Feshbach resonance. In other words we will take for 
the effective interaction the same as the one we had only two atoms present, namely Too = \w\'^ / {uj — Eq) as it 
results from Eq.(5). This is equivalent to retain only ladder diagrams for the short range potential. We note that 
this effective interaction is analogous to the one due to phonon exchange in standard superconductors, although there 
are differences. Actually we believe that this description should still be correct even if we take into account the 
perturbation due to the other atoms, provided it is not too strong. In this case we will still have resonant scattering, 
with possibly a shift in its energy position, a change in its width and its coupling strength. This can be accounted 
for by a modification of the effective parameters of the resonance. Moreover we will already find in our description a 
modification of the position and of the width. So the small perturbation due to the other atoms will possibly give a 
small quantitative change, but not a qualitative one. Naturally uj in our expression for Too is the energy for the center 
of mass of the two atoms at rest. If their total momentum is K and their total energy w we have to discount the 
energy associated with the center of mass motion and write roo(f^) = — Eq), with Q = w + 2jU — K'^/Am. In 

this last expression we have just taken into account that we will take in the following the origin of the single particle 
energy at the chemical potential /i, which produces a change for the energy of two atoms. 

We will now explore the simplest consequences of this interaction by ignoring fluctuation-like effects and analogous 
terms, and taking merely Foo as irreducible vertex. With this assumption we can write the Bethe-Salpeter equation 
for the full vertex V{lo,K), which is directly related to the scattering amplitude, as: 

r(a;, K) = roo(r!) + roo(f^)n(a>, K)T{uj, K) (9) 

where Yi{ijj,K) describes the propagation of two atoms and is given, in terms of the full thermal propagator of an 
atom G{uj, k), by: 

n{uj,K) = -TY, j ^GHa;-a;„,K-k)GK,k) (10) 

with LOn = (2n + 1)ttT being the Matsubara frequency. We have written Eq.(9) in a simple way by taking already 
into account that the full vertex F depends only, in our case, on the total energy uj and the total momentum K of the 
scattering atoms. This results from the fact that the irreducible vertex Foo(ri) has itself this property. We can more 
simply rewrite Eq.(9) as: 

V^^{VL)=T-^{uj,K) + I].{uj,K) (11) 

We will now eliminate the pole location Eq in Fqo, which is not an observable quantity, in favor of the physical energy 
Wo of the resonance, which is directly related to the scattering length as we have seen. This is done by writing Eq.(ll), 
at zero temperature T, for the case of two atoms in vacuum (implying /i = 0), at zero energy w = and momentum 
K = 0. In this case G(a;„,k) becomes the free propagator (icij„ — tk)~^ and the frequency summation in 11 becomes 
an easy integration. On the other hand for two atoms in vacuum T{oj,K) becomes ro(fi) = |t/;|^/(f2 — ujq), which is 
essentially the T-matrix. Hence in this particular situation Eq.(ll) reduces to: 

-E,IH^ = -.M'-j^,^^ (12) 

Actually this equation is just equivalent to u>q = Eq -\- ReJE'c- Subtracting Eq.(12) from Eq.(ll) we obtain: 

T-\,j,K) = V^\n)-n{u,K)- [ -^J- (13) 

Taken together the last two terms of the r.h.s. of Eq.(13) give an integral which converges for large values of k, while 
each term separately is divergent. However this divergence is not a real problem. It would not be present if we had 



3 



kept the k dependence of < kjfl^s > instead of making k = at the outset. Eq.(13) makes clear a general feature, 
namely the scattering amplitude depends not only on the total energy lu of the two particles, but also on their total 
momentum. This is due to the term n(w, K) which gives the effect of the other fermions on the scattering process. 
In other words Galilean invariance for this process is obviously lost because of the presence of the Fermi sea. It is 
clear physically that the existence of the Fermi sea wil be unimportant when K is very large. On the other hand we 
expect this feature will be most important for K = and we will consider only this situation in the following. 

Quite remarkably the modifications produced by the other fermions are already very important when the effect of 
the interactions is omitted. This corresponds to the modification of the scattering due to Pauli exclusion. We will 
restrict ourselves to this particular problem in the rest of the paper. In this case we have G{ujn, k) = (iiOn — Cfe + 
and the calculation of Il(oj,K) can be reduced to a single integration over the momentum k. It is convenient to 
use reduced units to display this result. We take /x as our energy scale and fco as a scale for wavevector defined by 
/i = kQ/2m (this is the Fermi wavevector at T = and not much different at low temperature). We introduce the 
reduced wavevector x = fc/fco, the reduced energy w = a;//i and reduced temperature f = T//x. ThenEq.(13) becomes: 

-mk,T(^)=-X-^^L '^^[l-^^3T3^tanh_] (14) 

where in the last term a) has to be understood with an infinitesimal positive imaginary part. Except for the factor 
7r/2fco this is just the inverse of the effective scattering amplitude for our problem. We have introduced the 
coupling constant A = 2fco|a|/7r, also related to the detuning wq by A = {2/it)W/ujq with W = 7fco being the energetic 
(half) width of the resonance line for a wavevector fco- We have used the reduced width W = {2/-K)W/iJb. If we consider 
the specific case of ^Li , the most heavily explored experimentally, the standard energy width of the Feshbach resonance 
is given [4] by 7/|abg| where af,g is the high field limit of the scattering length, of order of 100 nm. It is experimentally 
of order of 100 G, which translates into an energy of 10 mK. If we consider dense gases for which we have fco|a(,g| ^ 1, 
this gives ~ 10 mK. Since we will have at most experimentally /x ~ lO/xK, we see that W ^ 10^. Since we are 
interested in reduced energy Q of order 1, this makes the second term of the r.h.s. completely negligible and we omit 
it from now on. The same is likely to be true in most useful cases of Feshbach resonance. Note however that this 
term is necessary if we want to find in the lower complex plane the pole corresponding to the Feshbach resonance. 

Let us call I{u)) the integral in Eq.(14). The imaginary part of the r.h.s. is easily found to be Im/(u)) = 
— (7r/2)i?tanh(a}/4f) with R = [1 + (D/2)^/^, and is plotted in Fig. 1. At zero temperature with no Fermi sea, 
this would give us back the imaginary part in Eq.(8). We see that, in addition to w = — 2 corresponding to zero 
kinetic energy, this imaginary part is zero for w = 0, that is at the chemical potential. This is expected on general 
grounds since injecting particles at this energy does not perturb equilibrium and so does not lead to decay. More 
generally the tanh(w/4<) can be understood as the decrease of the scattering resulting from Pauli exclusion on the 
final state, together with the existence of reverse processes, both due to thermal occupation. 
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FIG. 1. Imaginary and real part of the integral I{uj) in Eq.(14), as a function of reduced energy u) = u)/n for various reduced 
temperature t = T/fj, indicated in the figure. 

Re/(a;) is plotted in Fig. 1 for various reduced temperatures. After adding the term 1/A and multiplying by the 

factor — 7r/2fco, we can consider the result as the inverse of an effective scattering length a~ff for two atoms. When A 
is of order of unity or larger, we see as naturally expected that the scale for this scattering length is the only one left 
in the problem, namely 1/ko- As seen in the figure a~j^{uj) has a strong energy dependence on the scale of the Fermi 
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energy Ep. This is in contrast with the case of two isolated atoms seen in Eq.(8), where in the same energy range the 
real part of — /^^ is a constant equal to the scattering length a^^ (except if luq is of order Ep that is extremely near the 
Feshbach resonance). This means that an essential simplification in the scattering properties effectively disappears. 
Indeed for ultracold atoms scattering is characterized by a single parameter, namely the scattering length a. Now 
because of the Fermi sea the energy scale Ep appears and the scattering amplitude gets a complex energy dependence, 
which depends also on temperature. An immediate consequence is that physically the Feshbach resonance is actually 
washed out by the Fermi sea. Indeed instead of having for all possible scattering atoms a diverging scattering length, 
and correspondingly a zero Re/~^, we have a Re/~^ which is of order 1/kp and depends on the energy of the two 
considered atoms (as well as their momentum as we have seen). This occurs as soon as A is not small. In particular 
nothing special occurs right at the Feshbach resonance when = 0. This provides a simple explanation to the 
experimental observations [2] that the resonance is not seen when the magnetic field is swept through its assumed 
location when the gas is dense enough to be in the degenerate regime. Naturally the inhomogeneity due to the varying 
trapping potential is an additional source of smearing since the energy scale /i for the scattering is space dependent. 

We look now more strictly for a resonance where the scattering amplitude diverges, which for two isolated atoms 
occurs at zero energy at the Feshbach resonance = 0. So we require that both the real and the imaginary part 
of r(w,0)~^ are zero in Eq.(14). It is easy to see that, for iD = —2, the real part of the integral in Eq.(14) is always 
positive so Re/^^ > for a < 0. So the only possible resonance occurs at the chemical potential uj = 0. In this case 
the condition that Re/~^ = in Eq.(14) coincides with the well-known condition for the BCS pairing instability. In 
particular, at T = 0, the logarithmic divergence of Re/~^ which occurs for a) = is a mark of this instability. It is 
known that this pairing instability is basically due to Pauli exclusion by the Fermi sea on low energy states, which 
produces a shift from a 3D situation to an effective 2D physics. We can understand qualitatively the strong energy 
dependence of the scattering amplitude we have found above as a manifestation of this 2D physics. Naturally, since 
we have not included interactions, the value of the critical temperature Tc we have for the pairing instability is just 
the standard one [3] , and it does not contain lower order fiuctuation effects [9] nor higher orders and self-energy effects 
considered in recent calculations [10]. Obviously our calculation is no longer strictly valid below Tc since we should 
deal with the superfluid state, which can be done as we have already mentionned. However this is not an important 
restriction for our present purpose since we want to stress the strong frequency dependence of Re/~^ on the scale 
of Ep, which will clearly remain at least rjualitatively valid in the superfluid state. In the same spirit we have not 
included interactions, but they will not change the basic Pauli exclusion physics. We can not expect interactions to 
remove the energy dependence of the scattering and our qualitative conclusion will remain unchanged. 

In conclusion we have presented a coherent framework which allows to deal with many-body effects in the presence of 
a Feshbach resonance. As a simple consequence we have shown that the mere result of Pauli exclusion, which results 
from Fermi statistics, induces a strong modification of the scattering properties. It is clear that this modification 
is a necessary ingredient in the physical understanding of these systems since Pauli exclusion can not be ignored. 
This modification results in a smearing of the Feshbach resonance and provides a natural explanation for recent 
experimental findings. 
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